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ABSTRACT 

This report presents the results of a theoretical and experimental 
investigation on the reduction of boundary layer thickness in low-density 
nozzles by wall cooling, wall suction, and a combination of these two. 
Potentially there is a twofold benefit in reducing the thickness of the 
nozzle boundary layer: (1) a possible increase in diffuser effectiveness, 

and (2) a possible reduction in the amount of boundary layer flow for a 
specified usable test section size, or an increase in the size of the 
usable test section for a given nozzle mass flow. The theoretical develop- 
ment starts with the proper integral relationship for a compressible 
laminar boundary layer. The normal velocity at the wall is allowed to 
be finite to include the effects of wall suction. Definitions of momentum 
and displacement thicknesses which account for transverse curvature are 
used. The results of Iglisch, who developed an exact solution for incom- 
pressible flat plate flow with suction, are used to estimate skin friction 
coefficients. The Prandtl number is assumed to be unity, and two-dimensional 
values of &*/ 0 are used. An exponential velocity profile which takes wall 
suction into account was used to estimate boundary layer height. 



The theoretical results were checked by an experiment in which a 
Mach number 9 - to - 11 porous nozzle was operated at unit Reynolds numbers 
in the range of 100/ inch to 600/ inch. Pitot pressure surveys were used to 
determine the exit Mach number and boundary layer thickness. Theoretical 
Mach number predictions are shown to agree with the experimental results 
to within 5%, and boundary layer height predictions to within 10%. 


Theoretical results are presented which show the effects of suction 
and wall cooling at several Reynolds numbers on nozzle diameter and uni- 
form core size for a given throat area and Mach number distribution. It 
is concluded that the use of suction and cooling may result in a larger 
test section size, but that the merits of a cooled porous wall in any specific 
case must be decided from an analysis of the complete wind tunnel system. 
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DEFINITION OF SYMBOLS 


Symbol 

a 

A 


Cf 

D 
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m 
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Pn 


P* 


Pr 

Re 




x, y, r, s 

7 

5 


Definition 
nozzle wall radius 
area 

local skin friction coefficient 

diameter of suction hole 

enthalpy per unit mass 

nozzle wall thickness 

mass flow rate 

pressure gradient parameter 

pressure 

stagnation pressure 

stagnation pressure behind normal shock 

Prandtl number 
Reynolds number 
temperature 

stagnation temperature 
velocity 

uniform suction velocity 
weight flow rate 

nozzle coordinates (see Figure 3) 
ratio of specific heats 
boundary layer thickness 
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DEFINITION OF SYMBOLS (Cont'd) 

Definition 

boundary layer displacement thickness 

boundary layer momentum thickness 

viscosity 

(v 0 /ui ) 2 Re! 

dens ity 

ratio of open area to total wall area 

ratio of effective open area to total wall area 

nozzle wall half-angle 

SUBSCRIPTS 

wall 

frees tr earn 
sonic condition 
adiabatic wall 
expansion core 
suction 
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SUMMARY 


This report presents the results of a theoretical and experimental 
investigation on the reduction of boundary layer thickness in low-density 
nozzles by wall cooling, wall suction, and a combination of these two. 
Potentially there is a twofold benefit in reducing the thickness of the 
nozzle boundary layer: (1) a possible increase in diffuser effectiveness, 

and (2) a possible reduction in the amount of boundary layer flow for a 
specified usable test section size, or an increase in the size of the 
usable test section for a given nozzle mass flow. The theoretical 
development starts with the proper integral relationship for a compres- 
sible laminar boundary layer. The normal velocity at the wall is allowed 
to be finite to include the effects of wall suction. Definitions of 
momentum and displacement thicknesses which account for transverse curva- 
ture are used. The results of Iglisch, who developed an exact solution 
for incompressible flat plate flow with suction, are used to estimate skin 
friction coefficients. The Prandtl number is assumed to be unity, and 
two-dimensional values of &*/ 9 are used. An exponential velocity profile 
which takes wall suction into account was used to estimate boundary layer 
height. 

The theoretical results were checked by an experiment in which a 
Mach number 9 - to - 11 porous nozzle was operated at unit Reynolds 
numbers in the range of 100/ inch to 600/ inch. Pitot pressure surveys 
were used to determine the exit Mach number and boundary layer thickness. 
Theoretical Mach number predictions are shown to agree with the experi- 
mental results to within 5%, and boundary layer height predictions to 
within 10%. 


Engineering Center, University of Southern California 


Theoretical results are presented which show the effects of suction 
and wall cooling at several Reynolds numbers on nozzle diameter and uni- 
form core size for a given throat area and Mach number distribution. It 
is concluded that the use of suction and cooling, may result in a larger 
test section size, but that the merits of a cooled porous wall in any 
specific case must be decided from an analysis of the complete wind 
tunnel system. 


I. INTRODUCTION 

One of the main problems in designing a hypersonic low-density 
wind tunnel is devising a pumping system that will handle the high 
volume flows at the tunnel exit. These high volume flows result 
primarily from the low pressure recovery associated with a diffuser 
operating with an entering flow field that is largely composed of 
boundary layer. A secondary cause of the high volume flow is the direct 
effect of the thick boundary layer in the nozzle. Since the boundary 
layer flow sometimes fills a large part of the nozzle, it is necessary 
in such instances to have a nozzle diameter many times that of the usable 
test section. This flow in the boundary layer must still be pumped, and 
it contributes significantly to the high volume flow. Thus, it appears 
that potentially there is a twofold benefit in reducing the thickness 
of the nozzle boundary layer: (1) a possible increase in diffuser 

effectiveness, and (2) a reduction in the amount of boundary layer flow 
that must be pumped for a specified test section size, or an increase 
in the size of the uniform test section for a given mass flow. 

This report presents the results of a theoretical and experimental 
study in the control o^ boundary layers in low-density wind tunnel 
nozzles. The work was carried out by the University of Southern 
California Engineering Center under National Aeronautics and Space 
Administration Contract NAS 8-5056. The NASA Marshall Space Flight 
Center, Huntsville, Alabama, instituted the contract in an attempt to 
develop an understanding of nozzle boundary layer control which would 
be useful in the design of nozzles for a low density test facility. 

The facility is to be of the cryopumped type (one-kilowatt refrigerator 
capacity), wherein the flow leaving the nozzle is condensed on an array 
of cold plates. The planned Mach number range is from 3 to 12. The 
maximum nozzle flow rate which can be continually accepted by a 1-KW 
cryopump is about 2.8 grams of nitrogen per second (this assumes that 
all the flow has been precooled to 100°K). Thus, although the boundary 
layer study herein described uses a general approach, the above considera- 
tions served more or less as guidelines, and the experimental portion 
of the study was performed in a small cryopumped tank. This tank closely 
simulated the conditions which should.be encountered in the full scale 
facility. 
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II. FACTORS AFFECTING NOZZLE BOUNDARY LAYER CONTROL 


For the purposes of this study, the "control" of the Boundary 
layer means reducing the height and preventing separation. For a 
fixed Reynolds number, the parameters having the greatest effect on 
the boundary layer height are (1) pressure gradient, (2) surface 
roughness, (3) ratio of surface temperature to frees tream temperature, 
and (4) surface suction or injection. These parameters are discussed 
in the following sections. 

A. Pressure Gradient 


In a nozzle, the pressure gradient is favorable (i.e., nega- 
tive), and increasing the pressure gradient will generally reduce the 
nozzle length and lower the boundary layer height in the test section. 
However, since the magnitude of this gradient depends upon the Mach 
number gradient, and large test section Mach number gradients are 
undesirable, there is a limit to the allowable pressure gradient. If 
the nozzle is to have a region of uniform parallel flow, there is a 
minimum nozzle length required; however, Reference 1 shows that the 
nozzle length can be significantly reduced if small Mach number gra- 
dients are acceptable. 

B. Surface Roughness 

Surface roughness has the potential of influencing the boun- 
dary layer height in two regimes. 

On the macroscopic scale, reducing the surface roughness will 
retard the transition to turbulent flow when the roughness elements are 
greater than some critical value. When the roughness elements are 
below this value, additional smoothing of the surface has no effect on 
the transition, and apparently the laminar boundary layer is unchanged 
also. With ordinary care the roughness of machined surfaces is below 
the critical value. 

On the microscopic scale, Reference 2 shows that surfaces are 
attainable where a significant fraction of the rebounding molecules 
*have their tangential momentum unchanged. In the limiting case with 
no loss of tangential momentum at the surface, the boundary layer could 
be eliminated; unfortunately, the techniques required to produce these 
surfaces are not suitable for wind tunnel construction, and the use of 
surfaces of this degree of smoothness does not seem feasible. 
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C. Temperature Ratio 


Reducing the wall temperature causes a reduction in boundary 
layer height. For example, Reference 3 shows that reducing the wall 
temperature from adiabatic wall temperature to freestream temperature 
reduces a two-dimensional laminar boundary layer by 50% at M = 7.0. 

The reduction is somewhat greater at the higher Mach numbers. Thus, a 
cooled nozzle wall appears to be a useful means of reducing the boundary 
layer height. This wall cooling also has structural advantages, since 
at high Mach numbers the adiabatic wall temperature is too high for most 
materials . 


D. Surface Suction and Injection 

There have been many papers published on the use of surface 
injection to reduce wall heat transfer, but the use of this process 
causes an increase in the boundary layer height. Reference 4 shows that 
the converse is true, that suction increases the heat transfer and reduces 
the boundary layer height. Stalder [5] describes the performance of a 
Mach number 4 porous nozzle developed at the University of California Low- 
Density Facility. This nozzle was constructed by bolting together a 
series of flat plates, machining the desired contour, and then removing 
every other plate to achieve the desired porosity. It was found that by 
varying the plate back pressure, the Mach number could be kept nearly 
constant over a wide Reynolds number range. 

To demonstrate the possibilities of suction and to present the 
overall trends, calculations were made on the reduction in boundary 
layer height that can be attained using suction and cold walls on a flat 
plate in a uniform stream. The computations used the equations of 
Reference 6 (summarized in Appendix A) and were made for a freestream 
Mach number of 10 with a unit Reynolds number of 600/ in and suction 
velocity of 2.26 % of the freestream velocity. These conditions are 
roughly those which one would obtain at a flow rate of 1 gram/sec in a 
12-inch diameter nozzle with a displacement thickness of 2.8 inches. 

The suction velocity is approximately that resulting from a 50%-open- 
area plate. The results of the calculations are presented in Figure 1. 

It is apparent that the cooling reduces the boundary layer height by a 
fixed fraction as the plate length is increased, whereas the suction 
causes an increasing fractional reduction as the plate length is increased. 
Thus the cooling produces another "similar" solution while the suction does 
not. The effect of the suction can be roughly visualized as imposing a 
uniform normal velocity on the flow over the plate. Since this imposed 
velocity is small compared to the freestream velocity, the boundary layer 
height is practically unaffected where the growth rate of the boundary 
layer is comparable to the freestream velocity, i.e., at the front of the 
plate. Aft on the plate, the boundary layer growth rate diminishes and 
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the imposed suction velocity is increasingly effective in reducing the 
height. This model is useful only for explaining the major trends and 
does not quantitatively represent the problem, however. 

For comparison, the boundary layer height has been calculated 
for a unit Reynolds number of 200/ inch, while holding the other factors 
fixed. The two cold-wall calculations are compared in Figure 2. Here 
it is seen that for a fixed longitudinal location, the fractional reduc- 
tion in boundary layer height is smaller for the lower unit Reynolds 
number case. This illustrates a possible difficulty in using this boun- 
dary layer control technique at extremely low unit Reynolds numbers. 

E. Combined Suction and Cooling in a Cryopumped Wind Tunnel 

In a conventional low-density wind tunnel, the feasibility of 
reducing the boundary layer height by means of suction is limited by the 
very large suction volume flow which must be pumped. In fact, an analysis 
by Enkenhus [7] demonstrated that the beneficial effect of suction on 
reducing the boundary layer height and increasing the test core size 
could be just as well accomplished by adding the suction pumps to the 
main test section flow pumps and increasing the nozzle size. 

There' is one pump, however, which can overcome many of the 
problems associated with high volume flows. This is the cryopump, which 
uses a highly cooled surface to condense the tunnel flow. In a typical 
installation using a helium refrigerator, extremely low pressures (of 
the order of 10“ 6 mm Hg or lower) are available if the tunnel gas is 
nitrogen or a combination of nitrogen and oxygen. These pressures are 
much lower than the static pressures of a typical low density facility, 
so it is possible to place the condenser plates in a chamber downstream 
of the test section and thus provide the low pressure region required 
to operate the tunnel. 

Rogers [8] has investigated the possibilities of boundary 
layer control in a small-scale wind tunnel of this type. It is shown 
in Reference 8 that the use of suction combined with highly cooled 
nozzle walls (liquid nitrogen was used in this case) results in a 
lower heat load to the cryopump, so that larger tunnel mass flows may 
be accepted. Thus, the use of wall suction and cooling in a cryopumped 
wind tunnel appears to be an especially promising means of reducing the 
height of the boundary layer and increasing the size of the test core. 

The effects of combined suction and cooling. are examined quantitatively 
in Section VI of this report. 
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Ill . THEORETICAL DEVELOPMENT 


A literature search was made for methods of computing the thickness 
of a laminar low-density boundary layer in an axisymmetric nozzle. The 
search yielded three approximate methods for boundary layer calculation, 
namely, those of Durand and Potter [9], Johnson [10], and Maslach and 
Sherman [11] . 


Of the above references, only the first takes into account transverse 
curvature effects in defining displacement thickness (5*) and momentum 
thickness (e) . The effect of transverse curvature on axisymmetric boun- 
dary layers is discussed by Probstein and Elliott [12]. This effect is 
important if the boundary layer thickness comprises a significant fraction 
of the nozzle radius, which is the case in a low density nozzle. The fol- 
lowing development therefore parallels the work of Durand and Potter [9] . 
Changes have been made where necessary to include the effects of wall 
suction. The coordinate system is given in Figure 3, where the expres- 
sion "expansion core" has been retained from Durand and Potter. This 
"core" is the one to which 6* is added to obtain the wall contour, but 
in the case of suction it does not determine the area ratio of the nozzle, 
since some of the flow has "expanded" through the walls. 


A. Boundary Layer Momentum Integral for Computation of 5* 

We begin with the continuity and momentum equations for axially 
symmetric flow. 

Continuity: 


JL 

5s 


(P“r> + ^ 


(pvr) = 0. 


( 1 ) 


Momentum: 




+ _d_ ( SuV & dr du 
ds dy \ Syy r <3y dy 


( 2 ) 


Multiplying equation (2) by r and integrating through the boundary 
layer from 0 to 6 gives 


6u , 

pur ^ dy 


dii , 

pvr^dy 



o o 


+ 



du dr , 

^ dy dy y * 


(3) 


6 


o 



Now if the continuity equation is integrated with suction 
allowed at the wall (v Q ^ 0) , it becomes 


pvr - p o v 0 a = 


y 


^ pur dy * 


( 4 ) 


Substituting equation (4) into equation (3), we have 


/ fr dy - / 57 [ j ' & 


pur dy 


dy+J^ P 0 v Q a dy 


o 

6 


* - / r & dy + / r |r( t ‘57) dy+ / ^ §7 §7 dy * (5) 


Integrating by parts the second term on the left hand side and 
the third term on the right hand side, and rearranging, we obtain 


/ Jr e“ 2r dy - “1 /-Jr / '<»■*» ( 37 ) 


o 

- / | 7 Po v o a d y- 


(6) 


If we substitute Euler's equation - dp/ds = PiU x for the 

pressure gradient, and integrate the last term of equation (6), we have 
(after multiplying by -1) 


du r 


f is pru2 dy + u i f ii pur dy + pi Ui §1 J r dy = a pqI 


(b uV 

w + u 


iPo v o a ■ 


( 7 ) 
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If the indicated differentiation in the first two terms is 
carried out, and the term 


6 

r 

J 


o 



dy 


is added to and subtracted from equation (7) 5 we get 



dy - 


pru 


Bs 


dy + 


dui 
PlUl Ss 


o 

I 


rdy 


= a |io 



+ u iPo v o a * 


o 


( 8 ) 


If we divide equation (8) by Pl u x 2 a, combine the second and 
third terms, and break the first term into two terms, we have 



-as— fi 

Pi U 1 \ 



dy + 


1 

Pi u i 2 a 



P u ■ 

Pi u i 




a 2 

^ pi u i a 



f?jn 
Pi «i 




Po Vo 
Pi U 1 


( 9 ) 


We now define the displacement and momentum thicknesses in 
terms of axially symmetric flow to account for the transverse curvature 
effect of a thick boundary layer. 


Mass Flow Defect: 


8 * 



Pl u-L dy = 


o 


5 


/ 


o 


2?tr (p x u x - pu) dy. 


(10) 
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Momentum Flow Defect: 


I 


o 


2 jar 


Pi Ui 2 dy 


8 


/ 


o 


2jtr p(uu x - u 2 ) dy. 

( 11 ) 


If it is now assumed that cos w = cos to Q = cos to ec t equations (10) and 
(11) can be written 


Si = 




( 12 ) 


where 


5i = 8* - 


8* 2 cos to 
2a 


(13) 


and 



(14) 


where 


01 s 0 


9 2 cos to 

2a 


(15) 


When equations (12) and (14) are substituted into equation (9), and 
dx cos to substituted for ds, we obtain the integral momentum equation 
for axisymmetric flow with suction: 


d0 

dx 


+ 


9l 


r (&-,/ e.i) 

L u i 


du 1 _d_ 

dx a p x u x 2 dx 


a 


Pi u l 



sec 


to + 


IQ. lx 

Ui T 0 


sec to 


where 


c f “ 


Au\ 

1/2 (pi ui 2 ) 


(16) 
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and 


(from assumption of no pressure 
change across boundary layer). 


Ii = m 
To Pi 


In place of the quantity 6^/ 0^, Durand and Potter recommend the 
substitution of (S*/e) 2 d. This is shown to be a very good approximation, 
and greatly reduces the complexity of the calculations. 

Equation (16) is a linear, first order differential equation, 
the solution of which can be given as 


01 » 


hJsl 

G(x) 




W(x) G(x) dx 


(17) 


where, with b>x/&i = (5*/e) 2<1 , 


G (x) = exp 


(5*/6) 2 d du 


Ui 


dx 


a Pi ui 


d . 
d£ (Pl 


a) 


dx (18) 


and 


W(x) 


J f sec co 


111 lx 

*1 T 0 


sec co« 


(19) 


If it is assuimed that the boundary layer calculations will be 
started at the throat with x = 0, M 1 = 1, and e x (o) = 0, then equation (17) 
becomes 

x 

01 “ GOO I W (x) G(x) dx. (20) 

o 

To solve equation (20) it is necessary to find methods of evaluating 
Cfj v 0 /ui, and (5*/9) 2 ( j* These items are discussed below. 


10 



Skin Friction Coefficient, Cj, and Form Factor, (6*/0 ) 2c i 


Lew and Fanucci [4] give an exact solution for uniform suction 
on a flat plate in compressible flow. They point out that when jj. is 
proportional to T, the equations have the same form in both the incom- 
.* press ible and compressible cases, and all the results for the incompres- 
sible case can be used. Appendix B shows the equivalence for the integral 
equations. 

Iglisch [13] has developed an exact solution for two-dimensional 
incompressible flow with s uction . Skin friction coefficient and (5*/ 0 ) 23 . 
are given as functions of si |£ nc The nondimens ional distance along the 
plate, | inc> , is defined as | in ^. = (vq/uj) 2 Re inc , where Re inc = poUjx/no. 
To obtain the corresponding compressible values in terms of. the freestream 
Reynolds number it should be observed that 


Re. = (Re ) 


me 


comp wall 


= (T3/T0) 2 Re 


comp 


since \j. ^ T and the static pressure is assumed constant through the 
boundary layer. Thus the results of Iglisch for Cf and (S*/0)2d* can 
be used directly as functions of 



instead of s[T~ m . 

*mc 


Figure 4 was taken from Table II in Iglisch. The value of (S*/ 0 ) 2 d^ nc 

obtained from the figure is used in the following equation (derived in 
Reference 14) to get (5*/0) 2 d : 

T 

(s*/e) 2 d = (&*/ e ^2d^ nc ^ " 1 ’ ( 21 ) 

If pressure gradient is to be taken into account, the results 
of Cohen and Reshotko [15] can be used. They give the following equa- 
tion for 5 */ 0 : 

6*/e = H t + 2 -= r ^- M x 2 (H t +1), ( 22 ) 

t-r z u. u r 

where is a function of wall to freestream enthalpy ratio and a pres- 

sure graSient parameter, n. Determination of the pressure gradient 
parameter is discussed in Reference 15. 
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The following equation for Cf , which includes the effect of 
pressure gradient, is given in Reference 15: 


c 

T jRe 1 ffi/To ) 2 -tj - • <23) 

(£ is a function of n and Tq/T^ 

The influence of pressure gradient was not accounted for in the present study 
Suction Velocity Ratio, Vo/u 7 

Since the velocity ratio is the dominant factor in reducing 
the boundary layer height by suction, it is necessary to relate this 
ratio to the pore or hole geometry used on the nozzle wall. In the 
theoretical analysis, it is assumed that a uniform normal velocity 
exists at the wall. In any physical case, it will be necessary to use 
an array of holes with discrete spacing. By making the hole size and 
spacing small enough, it should be possible to obtain a nearly uniform 
normal velocity. 

If the gradients in the boundary layer have a negligible effect 
on the flow through the holes, the holes may be considered as tubes with 
reservoir conditions equal to the local wall static pressure and local 
wall temperature. If the tube exit pressures are assumed to be negligible, 
the problem becomes one of determining the "choking" mass flow through 
tubes of varying L/D and Reynolds number. In the design of a nozzle, the 
length "L" will be the wall thickness which will be determined from a com- 
promise involving considerations of thermal conductivity, structural 
length, and the suction mass flow. It is desirable to have the walls 
thick to assure adequate strength and to minimize the number of cooling 
coils that must be used on the walls; however, at low Reynolds number 
it is necessary to minimize the wall thickness to allow high suction 
mass flows with the attendant high suction velocity ratios. 

There is one regime where simple, accurate estimates can be 
made of the mass flow through tubes. This is the free molecular flow 
regime, where it is assumed that the gas density is so low that the 
molecules collide with the walls many times before they collide with 
one another. In this limiting regime, the effects of L/D can be readily 
assessed using published graphs (Reference 16, for example). 
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For the particular case of a sharp-edged orifice, the con- 
tinuum mass flow has been computed theoretically and roughly checked 
experimentally, as shown in Reference 17. Here it is shown that with 
adequate pressure ratio (greater than 13) the ratio of actual mass flow 
to that existing with a uniformly choked opening is about 85%. Further, 
this reference presents the experimentally determined change in mass 
flow ratio as the Reynolds number is reduced to free molecular condi- 
tions. This data was limited to L/D = 1/40. Brown, et al. [18] give 
experimental results at higher values of L/D. Data from both Refer- 
ences 17 and 18 are shown in Figure 5, along with free molecular values 
of various L/D. The Reynolds number, Rej), is based on sonic flow. 

While the data in Figure 5 all have similar levels at small 
values of the parameter 64/Re D L/D, at the higher values of this param- 
eter the data show a definite trend with L/D. Thus at the higher values 
of 64/Re£ L/D, tubes with larger L/D ratios have lower mass flow ratios. 
The final value attained is the free molecular flow value, and the curve 
then ceases to change with 64/Re D L/D. 

Using Figure 5 with suitable fairing, m/m^. is known and v Q is 
given by the equation (negative sign for suction) 

*■*(%) 

* 

v ° = ■ Po Mo 


where 


m g is the suction mass flow rate, 

mo is the suction mass flow rate with uniform sonic 
flow through the holes, 

Po is average density at the wall, and 

M 0 is wall area 

for the nozzle Ax considered. Using standard equations for sonic mass 
flow (Reference 19, p. 82), the ratio v 0 /u! can be expressed as 

Vo/ux = - .0165 JJ -J%7 Tl (24) . 

where 0.is the ratio of open area to total wall area. 
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B. Determination of the Boundary Layer Thickness 


All of the discussion up to this point has concerned computa- 
tion of the displacement thickness, Mach number, and nozzle radius. 

Of equal importance is the determination of the actual thickness of 
the boundary layer, 8. It is this quantity which determines the size 
of the uniform test core. In this development, 8 is defined as the 
height where the velocity has reached either 99% or 99.9% of the free- 
stream velocity. The equations are set up so that either of these or 
any other desired velocity ratio can be specified. 

Equating Equations (12) and (13) gives 


8 * - 


8* 2 cos 03 
2a 



o 



_£U_\ 

Pi 


dy. 


(25) 


Substituting r = a - y cos co, and using y e for 8, where the 
subscript e is used to denote the edge of the boundary layer, we have 


/* 0 - ^ t C 1 - ^ 

o o 

(26) 

For ease of computation, equation (26) has been further simplified by 
breaking up the integrals as follows: 


5 * - 


5* 2 cos co _ 

2a y e 




J 


_£U_ 


Pi »1 


dy 


COS U) 

a 




- P u 
Pi ui 


y dy 


(27) 


Evaluation of the integrals in equation (27) will now be 
discussed. The development draws heavily on the work of Lew and 
Romano [6], which has been summarized in Appendix A. The first 
integral may be written 


X = 



- P u - 

Pi u i 


dy 



P u - d Z 5 
Pi ^l dt t 


dt. 


(28) 


14 



where the variable of integration has been changed from y to t (see 
Appendix A). It is shown in Appendix A that 


dZ = 

dt 



With this substitution, we have 


X = 



(29) 


Substituting for 5 t the expression 5 t = 2A/Re v Q (see Appendix A) and 
using the definitions of Re and v D given in Appendix A, equation (29) 
simplifies to 


x _ 2A u.n Tj. 
Vo Po T 0 




. 2 A MO 
V 0 Pi 



o 


o 


(30) 


To carry out the indicated integration, it is necessary to 
assume a velocity profile. In Reference 9, a velocity profile of the 

form u/u! - sin is used and good correlation with experimentally 


obtained values of 5 is shown. An exact solution for the suction case 
is given in Reference 4, but the computations are laborious. Lew and 
Romano suggest an exponential profile of the form u/u x = 1 - e“ (1 - tK) 
(see Appendix A). This profile showed good agreement with an exact solu- 
tion for the two-dimensional asymptotic suction case [4] at Mach number 
10 and Tq/T^ = 4. For the same conditions, it over-estimated the boun- 
dary layer thickness for the no-suction case (exact solution. Reference 20) 
by about 20%. Substituting the above velocity profile (equation (11) of 
Appendix A) and integrating, we get 


x = |?UiU [ T - 1 + K], (31) 

v o Pi e 

where terms involving e~ T e and e~ 2T e have been dropped since they are 
negligible. A is determined by Lew and Romano's solution, equation (13), 
Appendix A. 
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(32) 


| = \ (To/T ^ 2 


A(3A + 5) 
A + 1 


+ 5 in 


(A + 1) 
(2 A + 1) 


5 


where 


= (Vq/u-l) 


2 ,R.l. - U ..l, 


M-i 


All quantities needed to evaluate g are known, and A is determined by an 
iteration of equation (32). With A known, K may be evaluated, and r e 
determined by equation (11) of Appendix A for u/u x = .999, or whatever 
velocity ratio has been selected to define the edge of the boundary 
layer. 


The second integral in equation (27) may be written (following 
a development similar to that for the first integral): 


Y EE 



. R u . 

Pi u i 


y 


dy 


Zh PO 

v o Pi 



u j 

— y dr. 

Ui 


(33) 


The expression for y as a function of t is obtained by integrating 
equation (6) of Appendix A, which may be written in the form 


y(r) = 


2A no Tp 
v o Pi 


T 





1 


+ 




2c p To 


o 

Substituting the exponential velocity profile for u/u x 
the integration, one obtains 


dT. (34) 

and carrying out 


(t) _ 2A up Tp 

v o Pi T i ’ 


(35) 
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where 



t + 2e“ T 


2Ke" T (t + 1) 


1 

2 


e -2T 


(1 


2tK - K) 


- K 2 e 2T 




(36) 


+ 



1 



t + e 


-T 


-T 

e 


(t + 1) K - 1 4* K 



If we now replace the integral expressions in equation (27) by 
equations (31) and (33), and use equation (35) for y, we get 


5* 




S* cos 
2a 


w 


= y. 


2A H’Q 

v 0 Pi 


[T e - 1 + K] 



(37) 


In this equation &*, a, |a 0 > v Q , px> To/Tx, and cos w are all 
known quantities. The term y e is given by equation (35) which becomes, 
after terms involving e” T e are dropped (x e is always of the order of 
4.5 to 8), 


^ J- l .o la 
v 0 Pi Tx 




+ 



1 + 


2 i 

Ml To. 


(T, 


- 1 + K) . 


(38) 


In the above set of equations, A and all things determined by 
it (namely, K, x e , (u/u x ) e and Q) must be considered as a function of a 
flat plate x, which we could call x' . They cannot be evaluated at the 
particular nozzle x under consideration, since then equation (37) would 
be overspecified. It is noted that A is determined by equation (32) and 
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the x f dependence enters as the Reynolds number. Thus, to determine 6, 
we assume for the moment that the boundary layer is on a flat plate with 
a known Mach number, T 0 /T ± and 5*, and vary the x f until equations (37) 
and (38) are simultaneously satisfied. The term 5 is, of course, the 
value of y e for which the equations are satisfied; 

One problem occurred in the computation of & by the method 
outlined herein. For sufficiently high unit Reynolds lumbers in the 
suction case, no solution for 5 could be obtained, since for the &*, 

Mach number and other conditions given (which are those occurring in 
a nozzle, not a flat plate), the boundary layer height on a flat plate 
does not reach a high enough value at any x 1 to satisfy the equations. 
Suction boundary layers do not continue to grow indefinitely, but reach 
a maximum height at the point where the suction velocity balances the 
rate of growth. This limitation was encountered only once in the course 
of this study. The conditions for which no solution could be reached 
occurred in a 10° half-angle nozzle at a Mach number of 7.66, a unit 
Reynolds number of 5200/ inch, with a suction velocity ratio Vq/ u x of 
-.0305 (the minus sign indicates suction). The wall- to- frees tream 
temperature ratio , T 0 /T l5 was 1.88. The limitation thus apparently 
occurs only at relatively high unit Reynolds numbers and suction veloc- 
ity ratios (the open area in the above case was 607, of the wall area) 
even for highly cooled walls. For most cryopumped wind tunnels, it is 
considered that the limiting conditions would rarely be encountered. 


IV. METHODS OF CALCULATION 
A. Displacement Thickness, 5* 

We now proceed to discuss the methods of solving equation (20). 
Two basic approaches are possible, depending on what information is given. 
In both methods, it is necessary to integrate equation (20) numerically. 
Also, an iterative procedure is necessary in both methods. 

The first method attacks the problem of calculating the Mach 
number and uniform core radius with a given nozzle shape and porosity 
distribution. The procedure is of the iterative type because the Mach 
number depends on the magnitude of the suction, and in turn the suction 
mass flow depends on the Mach number distribution. It is necessary to 
assume hole diameter D as a function of x. It is considered that the 
hole diameters can be safely made to be 20% of the expected boundary 
layer thickness. The calculations are not sensitive to the hole size 
chosen, and thus conservative values should be used. The nozzle is 
divided into several increments of length. Ax. The calculation is not 
particularly sensitive to the number of increments chosen. For Mach 
numbers of the order of 12, twenty increments usually gives an accuracy 
consistent with the accuracy of the theory. The calculation is begun 
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at the first x and proceeds down the nozzle, 5* and Mach number being 
calculated at each x. It is necessary to compute 5 only at the last x. 
To begin the calculation, it is necessary to assume a value of the Mach 
number at the first station considered. In the initial portion of the 
nozzle, where the boundary layer is thin, a good first guess is the Mach 
‘number which would occur with no boundary layer and no suction. Further 
down the nozzle, it is better to add small Mach number increments to the 
Mach number at the previous station until agreement is obtained. 

The calculation is summarized as follows: For x^, assume a 

Mach number M x , where n refers to the nozzle location. Compute p l9 T x , 
Pi, u l9 from known stagnation conditions and isentropic equations. 


Tt f Tt V Z ' 1 ref + 

^ ref T ref ° r ^ ^ T ref/ \ T i + S / ^ef 

Note : If throughout most of the nozzle 

the static temperature is below 180°R, 
the linear law will give the better result 
for nitrogen. 

Re x = 

M-X 



Re, 


D 


R ^ 


_D_ 

H* 




A2(7-D 


= Reynolds number based 
on hole diameter and 
sonic conditions. 


where p* is given by jj* = p ref - 


ref 




p 


1 avg. 



= average pressure over the interval Ax. 


0*/0 = same as m/nu. and is read from Figure 5 at proper value 
of 64/ Re^ L/D. 
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Vo/ujl = -.0165 n/r/ 7 n/T 0 /Ti 


^aw = (^- + Mi 2 n/Pt^ = adiabatic wall temperature, 


I = (vo/u ^ 2 Re-L 


C £ (J! 

-r- •si Re-, is taken from Figure 4 at proper value of n/| ^ 

Z lo 

(5*/e)2d^ nc is taken from Figure 4 at proper value of v | ^ 


(S*/6)2d = (B*/6)2d. (To/Tj.) + ^ - 1 

inc l ■» 


r(&*/0) 


(5* /e ) 2d 1 


Ui 


2d 


r(s*/e) 2d i 


lx 

n 


Ui 


-5X 


n-1 


l avg 


A*i = (^i) x - (^i) x 

n n-1 


( a Pi u i 2 ) 


G(x) = 


x 

n 


( a Pi u i^") v _ 


x=0 


exp 


x 


n 


r(S*/0) 


Z&t-i 


2d 


ui 


avg. 




= ± + ±-L — 

' ' 2 COS CO U-L T 0 COS 


CO 


[G (x) W (x) ] 


[G(x) W(x) ] + [G(x) W (x) ] 

X X - 

n n-1 


avg. 


Note: When x n _^ = 0> take G(x) W(x) = 3 (Cf/2 cos co) , 

Cfa 

since Cf = — an( j G(x) is of the order of 1, 




** This equation results from forming the ratio of the flow per unit 
area at any nozzle station to the flow per unit area at the throat. 



where m gc is the mass flow in the expansion : core at the station con- 
sidered. Thus, 
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If M x assumed is within a specified tolerance of the Mach number given 
by A/A^, the calculation of 6* and Mach number is complete. Otherwise, 
a new value of M x is assumed and the calculation repeated. 

The second method begins with a known Mach number distribution 
and computes the required wall radius and uniform core size. It is 
necessary to specify either the porosity distribution (open area to total 
wall area, even though the wall area is not known) or the suction mass 
flow distribution, m s /m^. Since the Mach number is known, all Mach number 
dependent items can immediately be computed. If it is the porosity distri- 
bution (percent open area) which has been specified, the calculation pro- 
ceeds in exactly the same manner as the first method up to the G(x) 
calculation. A value for the wall radius is now assumed. Then 0!, Si, 
and m s /m^ are evaluated as before; r ec is obtained from the equation 


IT" (r ec /r *> : 




* 

11 - - 4 / 

m, 


The value of 6* is given by 


6 * = 




2 


r 

ec 


cos 0 ) 


&1 . 


The wall radius a may now be computed from a = r ec + 5* cos go. This new 
value is now used throughout the computation and the process repeated 
until the computed a agrees with the assumed a to within a prescribed 
tolerance. In all computations performed during this study, this pro- 
cedure always converged. 

If it is the suction weight flow ratios which are specified, a 
wall radius is assumed and Vq/uj computed from the equation 


RM T 0 

v ° /u i ■ - u, av 'pi ■ 

where An s and Mo are the incremental suction mass flow and wall area 
for the Ax considered. 0! and are computed as in the first method; 
r ec , 6' c , and a new a are computed as above, and the same iterative pro- 
cedure followed. After a is found, the required 0 is given by 
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0 , M s 

0V0 Pi rz 

0 1 v R 



Z 


^_A2(7-l) 

2 yj 


where 0 f /0 is the same as m/m^ and is known from Figure 5. 

B. Boundary Layer Thickness, 5 

With 5* and Mach number known from either of the above two 
methods the thickness of the boundary layer is computed by simultaneously 
satisfying equations (37) and (38). In these equations 6*, a, jj.o, v 0 . Pi, 
T i/ To, and cos oo are all known. The computation proceeds as follows: 

(1) Assume an x* . 

(2) Find A by an iteration of equation (32). A table of 
values of this function is presented in Table II. The 
assumed value of x f is used for x in equation (32). 

< 3 > K - - 1/2 (frfj) 

-T 

(4) x e is determined from (u/u x ) e ■ 1 - e (1 - t^K) . 

The value of (u/u!) e is that which defines the boundary 
layer thickness (usually taken as .999, but in this pro- 
cedure it is arbitrary). 

(5) A value for y e can now be computed from equation (38) and 
substituted into equation (37). The integral in the last 
term of equation (37) can be evaluated either numerically, 
graphically or analytically, xx/n^ is given by 

1 - e (1 - tK) and Q by equation (36). Because of the 
large number of terms in this expression, a graphical or 
numerical integration procedure is recommended. 

(6) If equation (37) is not satisfied with the assumed x 1 , a 
larger value is assumed and the process repeated; S is the 
value of y e for which equation (37) is satisfied. 
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In the case of no suction, (v Q = 0), the term 


2A up 
v o Pi 

is indeterminate since A = 0 when v Q = 0. An expansion of equation (32) 
for the case when A 0 shows that 

i 

2A [xp 
v o Pi 


can be replaced by 


1,2 x 1 



Pi 


V. EXPERIMENTAL INVESTIGATION AND COMPARISON WITH THEORY 
A. Experimental Arrangement 

An experimental program was devised to provide a check on the 
theoretical calculations.* A sketch of the experimental layout is pre- 
sented in Figure 6. During operation, the nitrogen used as the test gas 
is obtained by vaporizing liquid nitrogen. This gas passes through the 
1.5-K.W. heater to the stilling chamber and then into the nozzle. Part 
of the flow passes through the nozzle pores to cryopump coils surround- 
ing the nozzle, while the rest continues through the nozzle into the 
vane-type "diffuser” that serves as a precooler. This gas passes 
between the vanes to the remaining cryopump coils. The cryopump is 
cooled by gaseous helium from a 350-watt refrigerator. A liquid-nitrogen 
cooled shield reduces the ratiation load to the cryopump. A roughing 
pump is used for initial tank evacuation, and a diffusion pinup is used 
to remove any noncondensible gases. 

The nozzle schematic is presented in Figure 7. The nozzle is 
a 13° half-angle cone fabricated in two sections. The forward section 


* Experiments were conducted in the Hyperaltitude Facility of the 

Environmental Division of the U. S. Naval Missile Center, Point Mugu, 
California. 
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was machined from solid copper and includes the contraction section, 
the throat, and the supersonic section to a diameter of 2". The 9/32" 
diameter throat is a cylinder about one diameter long with the corners 
slightly rounded. 

The second section is the porous section that extends from the 
2" diameter to the 12" diameter. It was rolled from a 3/16" sheet of 
copper after the pores had been drilled. The hole pattern and the result- 
ing porosity are also given in Figure 7. The porosity distribution is 
essentially linear with axial distance. The hole radii were selected to 
keep the ratio of the hole radius to local boundary layer height equal 
to about .1; however, near the front of the porous section this ratio 
approached .6. 

The liquid-nitrogen cooling coils were soldered on this section 
with approximately 6" between coils. The coils did not cover any of the 
pores . 


After joining, the two sections were given the final machining 
so the nozzle coordinates are within 1% of a true 13° half-angle cone. 

B. Ins trumen ta t ion 


Test section pitot pressure was measured on an Alphatron 
(NRC Model 520). A calibration made prior to the test showed that in 
the range of most of the testing (lOOp - 300p), the instrument was as 
accurate as the scale could be read. This corresponds to from ± 5 % 
to ± 1.7%. Limited data taken below lOOp required corrections of less 
than 10% to achieve reading accuracy. 

Stagnation pressure was measured using a Bourdon gage at pres- 
sures above 20 mm Hg and a McLeod gage for lower pressures. The 
accuracy of this system was approximately ±2%. 

Tunnel static pressure was measured using a thermocouple gage 
designed for the low-micron range. No calibration was available, but 
the reading accuracy and the repeatability were approximately ± ,2p in 
the testing range. 

C. Experimental Results 

1. Porous Wall Nozzle 

Radial pitot probe surveys were made for a range of stagna- 
tion temperatures and pressures. In general a .312" O.D., .270" I.D. 
pitot tube was used. Since these dimensions are the same order of magni- 
tude as the mean free path of freestream gas, tests were made on a larger 
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probe (.875” O.D., .785" I.D.) and a smaller probe (.125" O.D., 

.106" I.D.) to determine the effect of probe size on the indicated 
pitot pressure. The results of these tests, shown in Figure 8, were 
used to correct the pitot probe surveys. Only the data taken with 
stagnation pressures less than .5 psia required corrections of 10% or 
more. 


These corrected surveys are presented in Figure 9. For 
these stagnation conditions, the nozzle developed Mach numbers ranging 
from 9.15 to 11.0 at freestream Reynolds numbers ranging from 100/ inch 
to 500/inch. At the lowest Reynolds number, the pitot pressure ratio 
is lower than would be expected. Despite the region of uniform pitot 
pressure (obtained with the .125" O.D. probe), it appears that the 
boundary layer has merged and the stream stagnation pressure is decreas- 
ing. Since the pitot pressure is ratioed to the upstream stagnation 
pressure, this would account for the low ratio. As will be shown in 
the next section, this is consistent with the theoretical calculations. 

At the higher Reynolds number the profiles show an 
increasing pitot pressure as the probe moves from the center. This 
is typical of low-density conical nozzles when they are operated at 
higher densities. The extent of this higher density regime can be 
estimated from the ratio of boundary layer height to test section 
radius (5/a) . This is illustrated in Figure 10, where the ratio of 
the pitot pressure at the edge of the core to the centerline pitot 
pressure has been plotted versus the ratio (5/a) for a series of conical 
nozzles. The data covers the Mach number range 4.5 - 11.0 and cone half- 
angles from 10° to 20°. From this plot it can be seen that as (5/a) is 
reduced below about .7, the pitot profile variations exceed 5%. Thus, 
to avoid transverse Mach number gradients, simple, conical nozzles must 
be operated at low densities in order that (5/a) exceeds .7. For higher 
density operation, it appears necessary to use a smaller cone angle or to 
contour the nozzle walls. 

2. Solid Wall Nozzle 

In order to obtain a direct comparison of a solid-wall 
nozzle and a porous-wall nozzle, the exterior of the nozzle was covered 
with aluminum foil. While this prevented any outflow, it did leave the 
interior nozzle wall rough due to the closed pores; however, it was 
felt that these closed pores would not present any greater disturbances 
than they had when used as a porous wall. 

The pitot profiles obtained from this solid-wall nozzle 
are presented in Figure 11. The corresponding porous-wall data are 
also presented. From a comparison of the two sets of data, it is clear 
that the main effect of the porous wall was to increase the Mach number 
while maintaining approximately the same boundary layer height. 
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D. Comparison of Theory and Experiment 


Theoretical calculations were made for the various test condi- 
tions, and the comparison between the theoretical and experimental Mach 
numbers and boundary layer heights is presented in Figures 12 and 13. 

For direct comparison the theoretical values of Reynolds number/ inch 
have been used as the abscissa for both sets of data. 

Figures 12 and 13 show that, with the exception of the lowest 
Reynolds number/ inch point, the theory predicts the Mach number to within 
5% and the boundary layer height to within 10%. Since these calculations 
were made using a 145°R nozzle wall, while in the experiment the throat 
block was at a higher temperature (about 285°R), calculations were also 
made to assess the effect of this temperature distribution. The results 
showed that the temperature distribution changed the Mach number by less 
than .5% and the boundary layer height by less than 2%. The boundary 
layer calculation on the lowest Reynolds number/ inch point showed that 
the boundary layers had merged. This agrees with the experimental 
results discussed in the previous section. 

Although the ratio of mean free path to pore radius varied from 
about .1 to 1.3, the suction mass flow was computed using free molecular 
flow values of m/m^ for all conditions. It is estimated that this 
simplification reduces the total suction mass flow by less than 10%. 

Since this calculated suction mass flow varied from 7% to 17% of the 
throat mass flow, the error is less than 2% of the throat mass flow. 

The experimental data with comparable stagnation conditions 
were also used to determine the suction mass flow by assuming that the 
difference between the exit Mach numbers for the solid wall and the 
porous wall was due to the suction mass flow reducing the exit mass 
flow. In this case the suction mass flow is determined from 


m out _ ^ (A/ A*.) solid wall 

m^ ’ (A/A*) porous wall 


The area ratios are determined from the exit Mach numbers. This pro- 
cedure* assumes that the effects of small changes of Reynolds number, 
velocity profile, and 5/a are negligible. A comparison of the results 
of these calculations with the theoretical values, presented in Figure 14, 
appears to justify the above assumptions. 
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In an attempt to learn more of the details of the flow coming 
from the pores, temperature and pitot pressures were measured at the 
exit of a 5/8” diameter pore. Due to the probable nonequilibrium con- 
dition of this exiting flow, and the unknown pitot probe errors, the 
results were only qualitative. The measured temperature was about 250°R 
and the pressure was about three times tunnel static pressure. Since 
the temperature and stagnation pressure in the boundary layer reach 
these values at distances of the order of one mean free path into the 
boundary layer, these values seem reasonable f The measured exit tempera- 
ture being above the nozzle wall temperature suggests that nozzle pores 
with small values of L/D cannot be used to precool the suction flow to 
wall temperature. For complete cooling, the pore length and diameter 
both must be sized relative to the tunnel wall mean free path if the 
tunnel boundary layer temperature rises significantly in distances of 
the order of the mean free path. 

To check the validity of the theory at higher Reynolds numbers, 
calculations were made on the cooled porous nozzle reported in Refer- 
ence 8. This 10° half-angle nozzle had a .104” diameter throat, a 
1.9” diameter exit, and a maximum porosity of 607 o . The comparison of 
the results of the calculations and the experiments is presented in 
Figure 15. The theoretical values of Reynolds number/ inch have been 
used as the abscissa. With the exceptions of the two highest Reynolds 
numbers where the nozzle flow had strong transverse gradients, the 
theory predicts the Mach number to within 3%. Similarly, with the 
exception of a few points, the boundary layer heights are in good agree- 
ment. The scatter in the boundary layer plots is believed to be due to 
the difficulty in defining the edge of the boundary layer on the experi- 
mental pitot profiles . 


VI. COOLING AND SUCTION EFFECTS ON NOZZLE AND TEST CORE 
SIZE FOR GIVEN MACH NUMBER DISTRIBUTIONS 

Machine calculations were carried out using the second theoretical 
method outlined in Section IV to determine the benefits of wall suction 
and cooling for a fixed Mach number. A Mach number distribution was 
assumed and the resulting boundary layer thickness, nozzle diameter, 
and uniform core thickness were calculated. 

A simple, rapid method of arriving at suitable Mach number distri- 
butions to use in the calculations was supplied by the first theoretical 
method, where a given nozzle shape is specified and a nozzle Mach number 
distribution is predicted. With a few trials using conical wall shapes, 
it was possible to predict Mach number distributions with exit Mach 
numbers close to 3, 6, 9, and 12, the range which was of interest in 
this study. Thus, although the distributions used were not aimed at a 
particular test section flow field, they do reflect smooth expansions 
which would closely approximate the flow in actual nozzles of the same 
test section Mach number, and are adequate for showing trends. 
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It was endeavored at all Mach numbers to pick a throat size such 
that the resulting uniform core was about six inches in radius. A 
nozzle length of four feet was used in all calculations, and in most 
cases resulted in reasonable nozzle wall angles. 

Figure 16 shows the nozzle wall radii, boundary layer thicknesses, 
and uniform core radii which were calculated for four Mach numbers. 

The Mach number distributions and throat sizes that were used are given 
in Table I. Nitrogen was assumed as the nozzle fluid. In all cases a 
wall temperature of 180°R was used, corresponding roughly to that for 
a wall cooled by liquid nitrogen. Results are given for the no-suction 
case and for the case where 22% of the nozzle flow is removed via wall 
suction. A unit Reynolds number range up to about 300 per inch was 
selected, since for the nozzle size chosen this corresponds approxi- 
mately to the maximum weight flow which can be accepted by a one- 
kilowatt cryopump. 

It is readily apparent from the figure that large reductions in 
the boundary layer thickness occur with 22% suction flow. For the 
Mach number 12 nozzle, S is reduced about four inches at a unit 
Reynolds numbqr of 160/ inch, corresponding to a flow rate of 2 grams/sec. 
The uniform core radii, however, are remarkabley close at all Reynolds 
numbers, so that the net effect of suction in all cases shown in 
Figure 16 is to reduce the nozzle diameter. 

What may be of more importance, however, is the beneficial effect 
of precooling in the suction case. It is difficult to exactly assess 
the benefits of the precooling on the nozzle size and core radii, since 
this is considerably influenced by the nature of the overall enthalpy 
removal system, in particular, whether a diffuser and precooler are 
used downstream of the nozzle and how effective they are. In this 
regard it is probable that the reduced boundary layer in the suction 
case would aid attempts to diffuse and precool the nozzle flow. A 
calculation has been made using simple assumptions about the system 
to provide an estimate of the possible increase in core size due to 
the use of suction and precooling. If it is assumed in the no-suction 
case that all of the enthalpy of the incoming flow is removed by the 
cryopump, and, for the suction case, 22% of the flow goes through the 
walls where this suction flow is all precooled to 180°R (liquid- 
nitrogen coolant), then for a stagnation temperature of 1460°R (Mach 12) 
we have 


suction 


1460°R 


w. 


no suction 


• 78 ^1460°^ + • 22 * h 180°R* 


235 cal/ gram 

.78(235 cal/gram) + .22(90 cal/gram) 


1.16. 
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Thus the test section area could be increased by 16% with a resulting 
7.5% increase in the uniform core size. Somewhat higher suction rates 
would raise this number, but it is doubtful that suction flows exceed- 
ing 30% of the throat flow are practical, especially if it is intended 
to cool the suction flow to the temperature of the wall. 

Figure 17 gives the results of calculations on the effect of various 
degrees of wall cooling on a, 5 and r uc for Mach number 12 and a unit 
Reynolds number of 800/ inch. This corresponds to a flow of ten grams 
per second for the nozzle size shown. The high Reynolds number was used 
in this plot to keep the wall mean free path to reasonable values in the 
adiabatic wall case. Even for this Reynolds number, however, the adia- 
batic wall mean free path is about 25% of the boundary layer thickness. 
Results are given for the no-suction case and for a suction case with 
0 max = .66. Again, large decreases in 5 are possible with the use of 
cold walls, but the uniform core radii are very close in the suction and 
no-suction cases for any wall temperature. A perhaps surprising feature 
of the plot is the increase in uniform core radii as the adiabatic wall 
temperature (no cooling) is approached. Of course, without cooling the 
nozzle radii get very large, and for the particular case shown here, this 
would result in a very large wall divergence even if suction were used. 
Moreover, it is uncertain how successfully such a thick boundary layer 
could be cooled after leaving the nozzle, and therefore a hot wall may 
still be inferior to a cold wall with the resulting increase in acceptable 
weight flow. Obviously, research is required regarding diffusion and 
cooling of flows with thick boundary layers before meaningful conclusions 
can be made concerning this point. 

Figure 18 shows, for a typical nozzle, the maximum ratios of open 
area to total wall area which are required to achieve various suction 
flow ratios. A linear suction rate has been used in all cases shown; 
i.e., starting at .4 feet do™ the nozzle, the suction flow is linearly 
increased and results in a total suction flow of either 11% or 22% of 
the throat mass flow. This also results in a nearly linear 0, the maxi- 
mum 0 in all cases occurring at the nozzle exit. Physically impossible 
0 1 s are indicated for the 22% suction ratio of Mach number 12. This case 
was used in Figure 16 for comparison purposes, but should not affect the 
validity of the indicated trends. 


VII. CONCLUDING REMARKS 

Based on a comparison of the results of theoretical and experi- 
mental investigations, it is concluded that this theory predicts the 
Mach numbers to about 5% and the boundary layer height to about 10% 
when the theory is applied to liquid-nitrogen cooled low density 
nozzles. The investigation covered the range of Reynolds numbers from 
150/ inch to 3000/ inch and the range of Mach numbers form 7.5 to 11. 
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When the theory was used to determine the effects of wall porosity 
and cooling on the radius of the uniform core, it was found that, for a 
fixed throat mass flow and fixed Mach number distribution, this radius 
was virtually unchanged. Since the cooling and wall suction both served 
to reduce the boundary layer height, this means that the cooling reduced 
S*/S enough to keep the boundary layer mass flow constant, while the wall 
suction removed mass flow at a rate that kept constant the sum of the 
suction mass flow and the boundary layer mass flow. While this showed 
there was no direct benefit from cooled porous walls, there were important 
secondary effects. One of these was the precooling of the suction flow, 
which would allow a higher mass flow (and larger nozzle) for a fixed 
cryopump size. Another advantage was that the smaller boundary layer 
meant a physically smaller and more manageable nozzle. Further, the 
smaller boundary layer might result in diffuser pressure recovery suf- 
ficient to allow conventional mechanical pumps to handle a significant 
fraction of the flow, permitting a larger total mass flow and a larger 
nozzle. Since these secondary advantages depend upon the overall system, 
it is concluded that the merits of a cooled porous wall must be decided 
from an analysis of the complete wind tunnel system. 
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MACH NUMBER 


X (feet) 


.2 

1.43 

.4 

1.64 

• 6 

1.80 

o 

• 

r— 1 

2.07 

1.2 

2.18 

1.4 

2.29 

1.6 

2.39 

1.8 

2.48 

2.0 

2.56 

2.2 

2.65 

2.4 

2.72 

2.6 

2.80 

K> 

• 

00 

2.87 

3.0 

2.94 

3.2 

3.01 

3.4 

3.07 

3.6 

3.14 

3.8 

3.20 

4.0 

3.25 

Throat 

Radii 

(feet) 

.284 


TABLE I 

TIONS FOR FOUR-FOOT NOZZLES 


MACH NUMBER 


2.21 

3.21 

3.77 

2.75 

4.13 

4.89 

3.16 

4.79 

5.70 

3.80 

5.75 

6.94 

4.05 

6.12 

7.45 

4.28 

6.45 

7.91 

4.49 

6.74 

8.34 

4.68 

7.01 

8.73 

4.85 

7.25 

9.10 

5.02 

7.48 

9.45 

5.17 

7.69 

9.78 

5.31 

7.88 

10.09 

5.44 

8.06 

10.39 

5.57 

8.22 

10.67 

5.69 

8.39 

10.95 

5.81 

8.54 

11.22 

5.92 

8.70 

11.47 

6.02 

8.82 

11.72 

6.12 

1 

8.95 

11.96 


.08 


.0351 


.0234 







TABLE XI 


A (3A + 5 


+ 


2A + 1 



-A 

f(A) 

.26 

.6816 

.27 

.7594 

.28 

.8446 

.29 

.9381 

.30 

1 . 0409 

.31 

1.1540 

.32 

1.2787 

.33 

1.4165 

.34 

1.5692 

.35 

1.7390 

.36 

1.9283 

.37 

2.1405 

.38 

2.3795 

.39 

2.6504 

.40 

2.9597 

.41 

3.3159 

.42 

3.7309 

.43 

4.2212 

.44 

4.8107 

.45 

5.5373 

.46 

6.4640 

.47 

7.7090 

.48 

9.5385 

.49 

12.8017 

.499 

24.1271 

.4999 

35.6081 
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APPENDIX A* 


COMPRESSIBLE LAMINAR BOUNDARY LAYER ON A FLAT PLATE WITH 
UNIFORM SUCTION AND WALL COOLING 

Consider the compressible flow over a flat plate having a uniform 
suction velocity v 0 and constant wall temperature T 0 . 


I I I I 1 I 

. Vo 


If the suction velocity is small compared to the free stream 
velocity u, the Prandtl boundary layer equations can be used to describe 
the flow. 

Momentum Equation 


Bu , 

pu a^ + 


p dy by 



Continuity Equation 


+ = o 

ax + ay U 


( 1 ) 


( 2 ) 


Energy Equation 

If the Prandtl number is assumed to be one, and the wall tempera- 
ture constant, the energy equation can be replaced by the following 
relationship: 


L = Iu 

To Tj. 



(3) 


* This appendix is a synopsis of the work of Lew and Romano, Reference 6. 
A prior work by Lew [23] concerning suction on insulated walls greatly 
facilitates the understanding of Reference 6. 
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Equation of State 


p = pRT. 


The continuity equation can be integrated to give 


pv - Po v O - “ 


s pudy - 


( 4 ) 


If this expression is substituted into the momentum equation, the 
Von Karman momentum integral may be obtained in the form 


d_ 

dx 


o 

J pCuu-L - u 2 ) dy - 


i 

Po v o u i " Ro 


(5) 


The x and y variables are now transformed as follows: Let 


dy = L 


and 


’l + 2- 

L u i 

(Ll i , _ 

u i 2 'N 

u 2 1 

\T 0 " 2 

cp T oj 

2 cp T q _ 


dx = L Re 

ds 



dt 


where L is the plate characteristic length and, 


Re = 


Po v 2cp T 0 L 
Po 


Equation (3) can be written in the form 



cp T 0 




( 6 ) 


(7) 


( 8 ) 
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where 


u = 


u 


and 


Ui = 


Ui 


y 2cp T 0 


y 2cp T 0 


Combining equations (6) and (8), we have 


dy 


L Pi 2 
2cp T 0 p / - 



L 


£l 

P 


II 

To 


dt. 


(9) 


Substituting equations (7) and (9) into (5), and rearranging, one can 
obtain the transformed momentum integral 



u 2 ) dt - Re u x v Q 



( 10 ) 


where S t is a measure of the boundary layer thickness and 


v 0 = 




\ J 2cp T 0 


Reference 6 suggests both an exponential and a quartic profile for sub- 
stitution into the integral relationship. However, the quartic profile 
leads to impossible profile shapes for large Reynolds numbers, and the 
exponential profile has been selected. The exponential velocity pro- 
file is 


— = 1 - e“ T (1 - tK) (11) 

u x 

where t = t/8 t and K must be determined by the boundary conditions. 

The boundary conditions which the velocity profile must satisfy 

are 


t — > 


t — > 0: 


u -> u l5 



d 2 u 
c)t 2 


-» 0 
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The fourth boundary condition was obtained by evaluating the 
momentum equation at the wall for the case when p is proportional to 
T. (If is not assumed proportional to T, the calculations become 
considerably more complex). 

Equation (11) satisfies the first three boundary conditions and 
the fourth gives 


K = 


1 (1 + 2A) 

2 (1 + A) 


where A has been substituted for the combination of variables 


Re v Q B 
2 


Insertion of the velocity profile into the momentum integral relation 
leads to the differential equation 


dA _ 2R e V ° 2 (4 A 2 + 4A + 1) (1 + A) ] 

ds u x 12 A 4 + 36 A 3 + 37 A* + 11AJ * 


( 12 ) 


This equation can be solved by integration to yield 


i 


2 Ot) 


A(3A + 5) 
A + 1 


+ 5 £n 


(A + 1) " 
(2A+ 1) 


where 



On Ui X 
Ml 


(13) 


and is a known quantity. 

Equation (13) thus is used to obtain A by a trial and error pro- 
cess. Quantities such as or 6, which can be expressed in terms of 
A, v 0 /u 1} T^/To, and the Reynolds number, can then be determined. 
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APPENDIX B 

SIMILARITY BETWEEN COMPRESSIBLE AND INCOMPRESSIBLE EQUATIONS 


Examination of the transformed momentum integral, equation (10) of 
Appendix A, has shown that it has the same form and boundary conditions 
as that for incompressible flow with suction. The Karman momentum inte- 
gral for flat plate flow with Prandtl number of one and having a suction 
velocity v Q at the wall is 


_d_ 

dx 


S(or oo) 


P (uu 3 


o 


U 2 ) 


dy - u lPo v 0 - po 



( 1 ) 


The boundary conditions which must be satisfied by the velocity 
profile which is used in equation (1) for incompressible flow are 


y - 0: 


u = 0, v 0 


du \ _ joq /c> fu\ 

& Jo PO \dPJo 


u 


u 


1 » 


Su 

-> 0 

oy 


3 2 u n 

dy 2 °* 


( 2 ) 


The second boundary condition is obtained by evaluating, at the wall, 
the momentum equation 


pu 



du _ & 
pV dy by 



The transformed compressible momentum integral is 


_d_ 

ds 



- u 2 ) 


dt - Re Ui 


- & 
V ° " ^St 


(3) 


where 
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t is a normal distance parameter 
S t is a measure of the boundary layer thickness 



ds = L Re dx 


u 


1 


Ul 

</2 c p T 0 


v 0 


v o 

^2c p T 0 


Po ^2cp T 0 L 

u = -- • Re = 

s/2c p T 0 

L = characteristic length of plate. 
Equation (3) may be written as: 


Po 



o 


(uu-L - u 2 ) d(Lt) - PqVo^i = Po 



(4) 


Letting a|t = Lt and the upper limit be co, we have 


Po 


_d_ 

dx 



o 


U 2 ) 


dijf 


Po u i v o 



(5) 


This equation has been transformed in y only. 


The boundary conditions to be satisfied are 



* This boundary condition was derived using the viscosity relation 

= T/Tq. Fortunately, for the wall temperatures of interest in this 
project (i % e. , those near the temperature of liquid nitrogen), this simple 
relation exactly satisfies the accurate Sutherland relation in the vicinity 
of the wall. The Sutherland law may be written 

(Footnote continued at bottom of next page) 
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Thus 5 the momentum integral and boundary conditions have the same 
forms in the transformed compressible as in the incompressible case. 

In both cases, the boundary layer form factor 5*/ 0 is a function of 
the nondimens ional distance along the plate, which has previously been 
defined. Because of the similarity of the y- transformed compressible 
equations and the incompressible equations, | evaluated at the wall is 
equivalent in the two cases. As shown on page 11 of this report, the 
incompressible and compressible £* s are then related by the equation 




me 


i± 

To 




JL. _ 

(J-o 


(T + s T 
•/Tq T + s Tq 


where s is a constant depending on the gas considered, and the subscript 
o refers to wall values. The term 


IT /Tq, + s\ 

V Tq V T + sj 


is of course 1 at the wall, and the rate of change of this term at the 
wall is 


d_ 




q + , 1 
+ s 



1 Tq_+_s 

2 T + s 


1 

n/Tq T 


IT To + s i /dt\ 

n / T 0 (T + s) a J o \^y) Q 



for s = To 

=» 100°K for nitrogen. 
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